Non linear sigma models and quantum spin systems 
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Microscopic models of quantum antiferromagnets are investigated on the basis of a mapping onto 
effective low energy hamiltonians. Lattice effects are carefully taken into account and their role is 
discussed. We show that the presence of an external staggered magnetic field modifies in a non 
trivial way the usual mapping onto the non linear sigma model, leading to the appearance of new 
terms, neglected in previous works. Our analysis is compared with Lanczos diagonalizations of S = f 
Heisenberg chains in a staggered field, confirming the validity of the single mode approximation for 
the evaluation of the dynamical structure factor. The results are relevant for the interpretation of 
experiments in quasi-one dimensional compounds. Microscopic realizations of SU(4) spin chains are 
also discussed in the framework of spin-orbital lattice systems. The low energy physics is shown to 
be described by sigma models with topological angle 8 in one dimension. This mapping strongly 
suggests that the one dimensional CP 3 model (with 8 — ty) undergoes a second order phase transition 
as a function of the coupling. 
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I. INTRODUCTION 



Quantum spin models in low dimensionality are cur- 
rently used to describe the magnetic properties of sev- 
eral materials including rare earths, organic compounds, 
transition metals and copper oxides. Available experi- 
mental data in these systems, especially magnetic reso- 
nance and neutron scattering, allow to establish the lo- 
cal structure as well as the most relevant features of the 
long range correlations. These studies generally show 
that magnetic materials can be accurately described by 
short range spin hamiltonians, at least in temperature 
regimes where the effects of disorder, the presence of 
spatial anisotropies or of dipolar interactions are negli- 
gible. The possibility to use simple spin hamiltonians to 
understand the physics of real systems has always been 
the drawing force for the development of more and more 
accurate methods to study the phase diagram of these 
models. Purely analytical techniques, like spin-waves ex- 
pansions or mean field theoriesa have shown to be quite 
accurate when magnetic ordering is present. Numerical 
methods, based oh series expansions! 2 ], Quantum Monte 
Carlo simulations^! or Density Matrix Renormalization 
Group (DMRG)q have been successfully applied also to 
finite temperatures, to frustrated models and to quasi 
one dimensional systems, where fluctuations play a key 
role in determining the physics of the model. A compli- 
mentary class of theoretical methods, which has been de- 
veloped for investigating the effects of quantum and ther- 
mal fluctuations in many body systems, is known as the 
semiclassical approach and includes different techniques 
which have been used to provide a physical interpretation 
to experimental and simulation dataa u. 

Semiclassical theories have been widely used in frame- 
work of quantum statistical mechanics since a long time. 
The mapping of quantum models, in particular spin sys- 
tems, onto classical effective statistical models has been 



extremely useful in understanding many relevant fea- 
tures of phase diagrams and the possible occurrence of 
quantum critical points. For instance, the one dimen- 
sional Ising model in a transverse magnetic field maps 
to the (anisotropic) two dimensional Ising model: its ex- 
act solution therefore provides a simple way to investi- 
gate the critical properties of the Ising universality classH. 
The quantum to classical mapping can be usually justi- 
fied on microscopic grounds only in the low energy and 
long wavelength limit, where the short range features of 
the original quantum model become irrelevant. There- 
fore, strictly speaking, the use of semiclassical methods 
in quantum statistical mechanics is restricted to regions 
characterized by a diverging correlation length and gap- 
less low energy excitations. However, the microscopic 
quantum hamiltonian is often assumed to be faithfully 
represented by its classical counterpart in a large portion 
of the phase diagram provided the coupling constants 
present in the classical lagrangian are suitably reMormal- 
ized due to quantum (short range) fluctuations!! This 
expectation has been beautifully confirmed by the study 
of Heisenberg models in low dimensionality. In partic- 
ular, exact solutions and numerical analysis of one di- 
mensional antiferromagnets found exponentially decay- 
ing correlations and gapped excitation spectrum in in- 
teger spin chains. Conversely, semi- integer spin chains 
turned out to be gapless with power-law correlations. 
This picture fully agrees with the conjecture by Haldane, 
based on the presence of a topological term in the semi- 
classical action of one dimensional spin chainsa. More- 
over, quantitative analysis of the three dimensional Non 
Linear Sigma Model (NLctM) found a phase diagram 
which compares favorably W.jtb,. experiments on two di- 
mensional antiferromagnetscO'El suggesting that semi- 
classical approaches can be directly used for the inter- 
pretation of experimental data. 

The semiclassical mapping for Heisenberg antiferro- 



magnets has been also employed to study the effects of 
an external magnetic field on the model. Of particular 
physical relevance is the role of a staggered field, which 
directly couples to the order parameter: Such a stag- 
gered field can be realized in certain quasi-one dimen- 
sional spin one compounds where, by lowering the tem- 
perature, rare earths magnetic ions undergo a Neel tran- 
sition to a three dimensional antiferromagnet which gen- 
erates_axi alternating magnetic field acting on the Ni 2+ 
chainsO. A similar mechanism has been also invoked for 
the interpretation of quasi one dimensional spin one-half 
chainscJ. However, the effective classical action used in 
the literature has not been explicitly derived from the 
quantum hamiltonian and we believe it should be re- 
considered. In order to extract quantitative information 
from this NLctM a-liSingle Mode Approximation" (SMA) 
is usually adopted^ but a detailed numerical study of its 
accuracy, in this case is still missing, even if analytical 
studieso suggest that it may be justified only for the 
transverse channel. In this paper, we present a micro- 
scopic derivation of the semiclassical action. We obtain 
an effective low energy theory different from the one com- 
monly adopted in the literature. This theory is then an- 
alyzed in the weak and strong field limit and the results 
are compared to Lanczos diagonalizations. 

Other magnetic materials where the use of semiclassi- 
cal methods may be.-suggestive are spin-orbital systems, 
like Ceo compoundscJ or transitional metal oxidesE3, in 
which orbital degeneracy is present. The special points 
where the model has an enlarged SU(4) symmetry are 
particularly importanJ—iu-prder to understand the phase 
diagram of the modeOllj and deserve a detailed analy- 
sis. We therefore apply the semiclassical mapping to the 
two microscopic models presenting this symmetry. We 
show that in one case the model maps straightforwardly 
to the CP 3 NLerM (with topological term in the one di- 
mensional limit) while in the other case, we have been 
able to carry out the mapping only in the special case 
of one dimension, where again we find the same semi- 
classical action at a different effective coupling. Interest- 
ingly, bpth-Jattice hamiltonians have been exactly solved 
in D=lE30 with very different results: the first model 
has broken symmetry and gapped excitations, while in 
the second case it is gapless and critical. We believe that 
this mapping provides an important clue to the under- 
standing the phase diagram of CP™ models, which are 
shown to undergo a phase transition as a function of the 
coupling constant. 



II. SEMICLASSICAL APPROACH 

In this Section we briefly review and generalize a 
method, proposed few years agecj, for the derivation 
of low energy effective actions in bipartite spin systems. 
The main advantage of this technique, with.respect to the 
original procedure developed by HaldaneEl, is the possi- 



bility to keep track of lattice effects and to make direct 
contact with other useful microscopic approaches, like 
spin wave theory. The method will be later applied to 
Heisenberg models in an external field and to spin-orbital 
systems. 

Be H the hamiltonian of a spin model on a bipartite 
lattice. For instance, the celebrated Heisenberg antifer- 
romagnet is described by 



w = 2J z2 Sfl ■ ' SflH 



(1) 



R£B 6 



where S are spin operatrs, the site index R runs over 
the sublattice labeled by B and 5 is a primitive vector 
on the lattice connecting nearest neighbor sites. In or- 
der to evaluate the partition function Z = Tr exp(— (3Ti) 
we adopt the usual coherent states formalism: we first 
split the interval (0, /3) in a large number N of Trotter 
(time) slices, then insert at each imaginary time t a res- 
olution of unity based on coherent states defined at each 
lattice site R. Here we follow the standard 0(3) nota- 
tion |fl(i2, t) >: In spin S models the coherent states are 
labeled by the unit vector $l(i?,£) and are characterized 
by the requirement that < il\ S \Sl >= SCI. These states 
may be explicitly obtained by a suitable rotation in spin 
space of the highest eigenvector of the S z operator. The 
partition function can be therefore written as: 



Vn(R,t)e~ s "" 





Seff = I dt[s 2 ^2Yl n ( R > *) • U ( R + 5 ' *) 

•'° ReB 8 

J2<n(R,t)\n(R,t)> 



(2) 



Up to this point, the underlying lattice structure is fully 
present in the functional form (|2|) and no approxima- 
tion has been introduced. However, the imaginary Wess- 
Zumino term 



< Sl\n >=iS 



(n x n) 

1 + fiz 



(3) 



prevents a simple classical interpretation of the effective 
action. In order to obtain a mapping onto a physically 
transparent classical statistical model, it is convenient 
to exploit the bipartite nature of the lattice by explic- 
itly tracing out the degrees of freedom defined on sub- 
lattice B. This procedure can be performed analytically 
because the variables we are integrating out are coupled 
only to the classical external fields tt(R, t) defined on the 
other sublattice (sublattice A). Therefore this step just 
requires the solution of the single site time dependent 
problem: 



-F\B(R,t)] 



TrU 



f —*KW 



K(R,t)=S-B(R,t) 
B(R,t) = sY J il ( R + S,t) 



(4) 



In terms of the local free energy functional F[B(R, t)], 
the effective action becomes: 

s eff = 22 F[B(R,t)} - dtJ2 < n(R,t)\n(R,t) > 

(5) 
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ReA 



Notice that now the effective action just depends on the 
field $1(A, t) on sublattice A. The solution of the problem 
defined in Eq. (|j) can be obtained within perturbation 
theory for slowly varying field configurations. This re- 
quirement is satisfied in the (semiclassical) large S limit 
and represents the only, low energy approximation we 
need to introduce in the evaluation of the effective ac- 
tion. We notice that the lattice structure is not involved 
in this semiclassical analysis and the long wavelength ap- 
proximation is actually unnecessary within our method. 
Explicitly, we obtain the following free energy functional: 

F[B(t)} = J dt |e (*) + r 00 (t) - j dt'T Q1 (t)T w (t')\ 
Fijit) = < Ui(t)\uj(t) > exp j j dt 1 [e,(0 - e,(t')]] (6) 

where 6i(t) is the i th instantaneous eigenvalue of K(i?,t) 
and \ui(t) > is the corresponding eigenstate. In order to 
obtain this expression we have assumed that the ground 
state of K(i?, t) \u (t) > is non degenerate at every time 
t. The terms shown in Eq. (g) are correct up to sec- 
ond order in time derivatives, as usual in semiclassical 
approximation. The long wavelength limit of the Barry 
phase term Too in A[B(£)] has been shown in Ref.tij to 
exactly compensate the Wess-Zumino contribution in the 
effective action (||). As a result, the low energy and long 
wavelength effective action, becomes a non-linear sigma 
model with, possibly, a topological 9 term coming from 
the residual space dependence of the Berry phase Too . 

Now we are in the position to apply this method to spe- 
cific lattice hamiltonians: The procedure we have just 
outlined requires the explicit solution of the instanta- 
neous eigenvalue problem defined by the on site hamil- 
tonian K(i?,t) (0), the evaluation of the terms ^ij(t) in 
Eq. (Q) and finally the substitution of the results into 
the form of the action (H) . 



III. HEISENBERG MODEL IN AN EXTERNAL 
FIELD 

It is now established that the ground state proper- 
ties and the low lying excitation spectrum of antifcr- 
romagnetic Heisenberg chains strongly depend on the 
value of the on site magnetic moment. Semi- integer spin 



chains have power law correlations and gapless spectrum 
while, for integer spin, correlations decay exponentially 
and a gap to all excitations is presented (A = 0.41048 
for S —r-X)- This different behavior was first found by 
HaldaneO via a semi-classical mapping onto the 0(3) non- 
linear cr-model. This conjecture was later confirmed by 
numerical studies based on Lanczos diagonalizationscj, 
DMRGEl and MC simulationsB 

Remarkably, these idealized models have also an 
experimental counterpart: neutron scattering exper- 
iment on quasi-one dimensional materials-, such as 
Ni(C 2 H 8 N 2 ) 2 N02(C104) (NENP), confirmed^ that the 
essential physics of these systems is well described by a 
simple S = 1 Heisenberg hamiltonian that couples neigh- 
boring spins antiferromagnetically and takes account of 
the single-ion anisotropy by including an on site term 
D (Sf) 2 . The predicted Haldane gap has been measured 
with great accuracy in these experiments showing good 
agreement with theoretical and numericalEOt£l predic- 
tions. 

The effects of a magnetic external perturbation on an 
antiferromagnetic chain, frequently gives rise to unex- 
pected interesting phenomena. This is the case of Cu 
benzoate, a quasi one dimensional S — 1/2 antiferromag- 
net displaying a gapped excitation spectrum in an applied 
uniform field. Oshikawa and AfflecktLj interpreted the ex- 
perimental findings on the basis of an Heisenberg hamil- 
tonian where spins are coupled to a weak effective stag- 
gered field. This microscopic model gives a field depen- 
dence for the gap which agrees with experimental data. 

More recently, interest has been focussed on the study 
of the effects of external fields on S — 1 systems. In the 
case of a uniform magnetic field, the lowest triplet excita- 
tion states are split into a transverse and a longitudinal 
mode and the gap closes at a critical field H c , where Bose 
condensation of magnons takes place23. The synthesis of 
compounds of the form R 2 BaNiOs (where R stands for 
a magnetic rare earth) allowed to study the effects of 
a staggered magnetic field on the quasi-one dimensional 
chain of spin one Ni 2+ ions. The magnetic moment of 
the Ni 2+ couples with the R 3+ ions that are ordered an- 
tiferromagnetically below a Neel temperature Tn (typ- 
ically 16 A" < T^v ^ 80 A). This three dimensional an- 
tiferromagnetic matrix generates an effective staggered 
magnetic field on the Ni 2+ chains whose intensity can be 
tuned by varying the temperature below Tjy. In this way, 
experiments have been able to investigate the effects of a 
staggered field on the Haldane g^Pri the staggered mag- 
netization and the susceptibilityt!l'El. 

Stimulated by these experiments, few analytical and 
numerical studies attempted a theoretical analysis of 
spin chaiua-iii staggered fields by use of semiclassical 
mappingsE3E-3 and DMRGl3. While qualitative agree- 
ment can be easily attained, some discrepancy still re- 
mains between the NLerM approach and DMRG findings, 
noticeably on the form of correlation functions. This cir- 
cumstance is rather surprising in light of the very nice 
agreement between the NLcM predictions and numeri- 



cal data for the pure Heisenberg chairc3. 

The strong field limit is particularly simple because, for 
every S, the ground state can be accurately described by 
a Neel state with gaussian transverse fluctuations. For 
such a problem, spin wave theory (SWT) can be applied 
also in one dimension giving a transverse spectrum of the 
form: 



e k = SJ(\K\/S 



2D)2_4 7 , 2 



(7) 



where H is the staggered field, D is the space dimen- 
sionality and 7ft = J2i=i cos ^i- To leading order, the 
transverse dynamical correlations in imaginary time pre- 
dicted by SWT have single mode character: 



S±(k,u>) = S 



2 \H\/S + 2(D- lk ) 



(8) 



The longitudinal correlations may be expressed as convo- 
lutions of the transverse ones, implying that no longitudi- 
nal branch of elementary excitations is present. Accord- 
ing to the SWT approach, in the strong field limit, the 
longitudinal gap saturates at twice the transverse one. 

Much more subtle is the weak field case, where quan- 
tum fluctuations strongly contrast the onset of a magnet- 
ically ordered state. In order to understand this limit, we 
re-examine the derivation of the NLerM for a spin S chain 
in an external staggered field on the basis of the method 
sketched in Section II. 

The microscopic hamiltonian we consider here is just 
the antiferromagnetic Heisenberg model in an external 
field H that can be taken either uniform or staggered: 

n = E S WE S ^ +H ] ±H -E S « (9) 
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perturbative treatment is justified only in the low-energy 
limit which can be physically accessed in gapless systems. 
In one dimension such a requirement is satisfied in half- 
integer spin chains while, for integer spins, it holds only 
in the large 5* limit. In fact, a perturbative renormal- 
ization group analysis, ip-pero external field, predicts the 
exponential dependences A ~ exp — (nS). This shows 
that the existence of the gap is a purely quantum effect: 
it vanishes on approaching the classical limit (S — * oo) 
where it is still justified to derive the effective action per- 
turbatively also for integer spin systems. Interestingly, 
in the strong field limit, the one body problem can be 
again easily solved by considering small oscillations of 
the vector f2 about the direction of the magnetic field. 
The resulting effective action, to quadratic order in the 
amplitude of the oscillations, exactly reproduces all the 
lowest order results of SWT, including lattice effects. 

Specializing Eq.mh to the form of the effective field 
( |ll| ) we obtain 



F[B(R,t)] 



where the conventions are the same as in ([[]), and upper 
(lower) sign refers to uniform (staggered) applied field. 
Following the derivation of the previous Section, we fac- 
torize the problem in the two sublattices (A and B) and 
write the effective action as: 

s e ff = - [ *V ( < to(R,t)\£i(R,t) >TSn-n(R,t) 

Jo ReA I 

+ ^F\B(R,t)] (10) 



ReB 



where now 



B(R,t) =S^2n(R + S,t)+H. 



(11) 



Therefore, we formally reduce to the same one-body 
problem as stated in Eq. (|J). This can be solved per- 
turbatively in the case of zero temperature (/3 — ► oo), 
and slowly varying effective fields B. The resulting func- 
tional F[B] is then correct up to second order in space 
and time derivatives. In this context we stress that this 







dtiS 



|rh(i?,t)| 2 
2\B(R,t)\ 



S|B(i2,t)| + roo[m] 

(12) 



where m = Mr. The first term comes from the time inte- 

|X5| 

gration of Eq. (||) and the second one is the ground state 
eigenvalue eo (t) . Now we perform the continuum limit of 
the expression (112) assuming that the relevant configu- 
rations Q(R,t) are slowly varying functions of space on 
the scale set by the lattice spacing a. To lowest order in 
spatial fluctuations we have B(R, t) = Sqfl(Ro,t) + H 
where Rq — R — x is a reference site belonging to sub- 
lattice A (x is the primitive vector pointing in the x di- 
rection) and q = 2D is the number of nearest neighbors 
of a hypercubic lattice in dimension D. In the weak field 
limit, we just need to keep terms up to second order in 

6B(R,t) =H + SX, y [n( J R+<$,*)-n(i*o,t)]. By expand- 
ing m(R, t) we obtain: 



5m(R, t) = m(R, t) - fl(R , t) 

SB(R,t) n(R 0l i) ■ 8B(R,t) 



Sq 



Sq 



+ 0(SB 2 ) (13) 



This leads to an approximation of the Berry phase that, 
to lowest order cancels the Wess-Zumino term in the ef- 
fective action leaving a residual contribution: 



r 00 H- r oo["ORo,t)] 

~ iS5m(R,t) ■ (m(R,t) x m(R,t)) 

~~H-(n(Ro,t)xn(Ro,t)) 

q 



(14) 



In one dimension, the usual topological term, arising from 
the spatial derivative of $1 present in SB, also appears 
besides the contributions shown in Eq. (pT|). By taking 
the continuum limit, we finally get the form of an effective 
action in which one half of the degree of freedom have 
been integrated out: 



S ef f= J dt J^J Yq (" ( *' t] + % H X ° (i? ' t] 



a 2 S 2 \V R Sl\ 2 - (1 =F 1)5H • ft(R,t) 



■2mSQ (15) 



The topological charge Q is non trivial only in D = 1 
where: 

Q = — dxdtft-(d x ftxfl) (16) 

As a result, we obtain an effective NLerM describing spin 
S chains in (weak) uniform (— ) or staggered (+) field. 
The microscopic derivation allows to obtain explicit ex- 
pressions for the bare spin wave velocity c = 2Sa\D and 
stiffness p s = S 2 a 2 ~ D which coincide with those already 
known at H = 0. While this.derivation reproduces known 
results for a uniform fieldtEl, it differs from the effective 
action usually quoted in the literature, where the exter- 
nal field only couples to ft via the Zeeman term. It is 
instructive to give a simple interpretation to the formal 
result we have obtained: a staggered field H can be writ- 
ten as the sum of a uniform field of the same strength 
minus a field twice as strong acting only on one sublat- 
tice (say the sublattice A). In this way, by tracing out 
the degrees of freedom living on sublattice B we obtain 
the same NLerM action appropriate for a uniform field 
with the addition of a Zeeman term — 2H • ft: This is 
exactly what we formally found in Eq. (|l5|). 

In order to investigate the low energy spectrum of the 
NLerM previously obtained (|15J), we resort to a simple 
single mode approximation: the constraint ft 2 = 1 may 
be lifted through the introduction of a Lagrange multi- 
plier A(H) and a linear shift of the field ft. In this way, 
the dynamical correlation functions in imaginary time 
acquire a Lorenzian form: 



m « cgS 2 /a 

°fc," fc,— uu 2 + 2cg\ + c 2 k 2 



(17) 



< S±S* _,.. > 



2cgS 2 /a 



k.uj —k 



(lu 2 - H 2 ) =f 2i\H\uj + 2cgX + c 2 k 2 



where A may be determined by a saddle point equation, 



or by fitting numerical data and g 



ca 



D-l 



I p s . Here the 



z axis identifies the direction of the external field. The 
poles of the correlation functions directly give the disper- 
sion relation of the model. Recalling that the NLerM de- 
scribes the spin degrees of freedom on a single sublattice, 
the wave vector k is defined modulo w. From expressions 
( |l7| ) it is apparent that the correlation functions have 
single mode character with different dispersion relations 
in the transverse and longitudinal channel: In particular, 
at each k (modulo ir) the transverse excitation splits into 
two different branches centered around the longitudinal 
dispersion: 



A L (fc) = v/A(H) 2 +c 2 fc 2 
A T (fc) = A L (fc)±|H| 



This form of the energy spectrum should apply to arbi- 
trary spin S at low energy and weak staggered fields. 
It definitely differs from the predictions of the usual 
semiclassical treatments which lead, within the same 
single mode approximation, to spectra of the form 
y/A(H) 2 + c 2 k 2 in both channels. Interestingly, also the 
exact solution of the ID S = \J2~XY model in a stag- 
gered field along the z axis shows£j a similar excitation 
spectrum, providing some support to our semiclassical 
analysis. The isotropic S = 1/2 Heisenberg model in a 
staggered field has been also analyzed-bw-,bosonization 
techniques and conformal field theoryt3E3. The stag- 
gered field has been shown to open a gap both in the 
transverse and in the longitudinal channel. This massive 
triplet is degenerate to leading order in the external field: 
Ay = Ai oc H 2 / z . It would be interesting to study the 
spliting of the gap by including subleading terms in order 
to compare conformal field theory results with the NLerM 
approach. 

In order to ascertain the validity of the single mode 
approximation in one dimensional models, we performed 
Lanczos diagonalizations in S — 1 chains. In particular 
we calculated the excitation spectrum by selecting, for 
each k, the excitation energy with the highest weight in 
the Lehmann representation of the dynamical correlation 
function. 
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FIG. 1. Longitudinal (above) and transverse (below) spec- 
tral weights for H = 0.1 (full dots) and H = 0.01 (open dots) 
staggered field. Lanczos diagonalizations in a N=18 lattice 
chain. In both cases, SMA is valid near k = n. 

In Fig. [y we show the largest normalized spectral 
weight Z(k) as a function of the momentum k of the exci- 
tation, showing that its value is always very close to one 
for k ~ 7r: this implies that the sum rule Q^ n Z n (k) = 1) 
is almost exhausted by a single excitation both in the 
transverse and in the longitudinal channel, thereby sup- 
porting the single mode approximation usually adopted. 
Notice that the excitation with largest weight does not 
always have the lowest energy, at fixed k. Lanczos diago- 
nalizations also show that the matrix element giving the 
spectral weight of the longitudinal excitation at k — n de- 



creases very quickly when the staggered field is switched 
on: It is lowered by a factor 2.5 when the field reaches the 
value H — 0.05. This result should be compared with the 
modest decrease of the transverse spectral weight, which 
reduces just by 20% in the same range. Both findings 
agree-,with the reported behavior of neutron scattering 
dataB. 

In Fig. g the energy spectra are shown for different 
field strength. While at low H the spectrum is markedly 
asymmetrical around k — tt/2, symmetry is restored at 
larger fields where it closely approaches the form ffl) pre- 
dicted by SWT. Data comes from Lanczos diagonaliza- 
tions performed on chains with 12, 14, 16 and 18 lattice 
sites. 
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FIG. 2. Transverse energy gap for two different field val- 
ues: H — 0.01 (open dots) and H — 1.0 (full dots). Lanczos 
diagonalizations in N=12,14,16,18 chains. In the high-field 
limit, the spectrum is well described by SWT (solid line). 

A specific feature of our results is the splitting of two 
branches in the transverse excitation spectrum when a 
weak staggered field is applied (|18|). Fig. |3| shows Lanc- 
zos diagonalizations data on several S = 1 chains which 
provide a numerical confirmation of the semiclassical pre- 
dictions. For every wavevector k we selected the longitu- 
dinal and the transverse excitation with the largest spec- 
tral weight Z. When the continuum limit is appropriate, 
i.e. at sufficiently small k (modulo n), two distinct ex- 
citation branches clearly appear, differing by ±|H| from 
the longitudinal excitation, in agreement with the NLctM 
analysis. The noticeable deviations around k — ir/2 are 
clearly due to lattice effects which are not correctly re- 
produced in the continuum limit. Notice that finite size 
corrections do not seem to affect the overall structure of 
the excitation spectrum of the model. 
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FIG. 3. Difference between the longitudinal and the trans- 
verse excitation spectrum for H = 0.01 and several chain 
length (N=12,14,16,18). The dashed lines correspond to the 
theoretical prediction Az, — At = dcH. 

In Rcf£3 the DMRG technique has been applied to 
investigate the lowest excitations in the longitudinal and 
transverse channels at k = n. The comparison of our 
results with DMRG is particularly instructive because it 
clearly shows the regimes where the use of the effective 
action approach is justified. 
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FIG. 4. DMRG dataO: we plot the difference between 
the longitudinal gap and the lowest transverse excitation at 
k = n. The data are compatible with our theoretical predic- 
tions in both the low and strong field limits. 

In fig. [| we plot the quantity 7 = (Al — At)/H 
as a function of H. The numerical results identify two 
different behaviors: at strong fields (H > 0.5) 7 satu- 
rates at 7 = 2, in agreement with SWT, while in the low 
field limit (H < 0.1) the splitting between the longitudi- 
nal and transverse gap quickly increases. In particular, 
the value 7 = 1 predicted by our semiclassical approach 



is compatible with the numerical data in the H — > 
limit, although extended calculations in lower fields are 
required in order to validate our analysis. 



IV. SPIN-ORBITAL MODELS 

Spin orbital models have recently attracted consider- 
able interest in the attempt to explain the unusual mag- 
netic properties of a class of quasi one dimensional mate* 
rials, which includes Cgo compounds (eg. TDAEj-C6o)lL3 
and few metal oxides (eg. Na 2 Ti 2 Sb0 2 , NaV 2 5 )ti. The 
physical properties of these Mott insulators are largely 
determined by the coupling between orbital and spin 
degrees of freedom which may be dominated either by 
Hund's rule or by dynamical Jahn- Teller effect. Possi- 
ble realisations in highet-dimensions are also found in 
fulleridescH or in LiNi0 2 H3. The low energy physics of 
these systems may be described by keeping only spin and 
orbital degrees of freedom. If the orbital degeneracy is 
twofold, like in the previous examples, the low energy 
model can be written in terms of two sets of spin- 1/2 op- 
erators per lattice site representing respectively spin (S) 
and orbital (T) degrees of freedom. Usually, spin isotropy 
in this effective low energy hamiltonian is retained only 
in the physical spin variable S while terms which break 
rotational invariance in the pseudo-spins T are generally 
allowed. However, the fully isotropic hamiltonian: 



(19) 



<i,3> 



<hJ> 



has been the subject of several studies, particularly in 
the two special cases K = ±4J. The K = +4J hamil- 
tonian can be written in terms of permutators on each 
lattice site and may be relevant for TDAE-Cgo while the 
K = —4 J naturally arises as the strong coupling limit of 
a microscopic hamiltonian appropriate when dynamical 
Jahn- Teller effect prevails. These two particular mod- 
els have several remarkable properties: besides an ob- 
vious SU(2)xSU(2) symmetry, thejt are both invariant 
by a larger SU(4) symmetry grouptLZI. The 15 generators 
of the symmetry group include the total spin and pseu- 
dospin operators: J2rSr an d J2rTr aim the further 9 
operators Y^r(^) R SrTr where the ± sign corresponds 
to the two models K = zLAJ. Note that, in the — case, 
the SU(4) generators do not commute with the trans- 
lations by one lattice spacing although the hamiltonian 
does not break any symmetry of the lattice. Remarkably, 
this model is also non frustrated: In a valence bond basis 
the ground state can be shown to have positive semi- 
definite weights. This feature allows to perform-very ac- 
curate Monte Carlo simulations on this systemEH. 

Both r&adeJs can be exactly solved in one 
dimensionsE3ty with very different physical properties: 
the K = +4J model (Sutherland model) is gapless, with 



power law spin correlations whose leading asymptotic be- 
havior has a x~ 3 / 2 decay and is characterized fw-,oscil- 
lations with period equal to four lattice spacingc3. In- 
stead, when K = —AJ (Valence Bond model) the system 
spontaneously dimerizes, the energy spectrum is gapped 
and correlations decay exponentiallyoO. Several rele- 
vant features of the ground states of these hamiltonians 
have been argued to be applicable to wider regions Jn 
parameter space, also outside the special SU(4) pointsLa. 
Here we will derive the effective low energy lagrangian 
for both models in order to understand how such a dif- 
ferent physical behavior in one dimension may originate 
and to shed light on the phase diagram of the two mod- 
els in higher dimensions. In fact, analytical and numeri- 
cal studies in two dimensions have suggested that SU(4) 
symmetry is not spontaneously broken .in, the ground 
state and a spin liquid phase may emergeL_fE3. The gen- 
eral method developed in Section II can be straightfor- 
wardly applied also to this class of hamiltonians defined 
on bipartite lattice because the interaction just couples 
nearest neighbor sites. Let us discuss the two cases sep- 
arately. 



A. Valence Bond model (K = — 4J) 

The model has four orthogonal states per site in 
the lattice, which correspond to the four possibilities 
(±^,±5) for the z-components of the spin and pseu- 
dospin variables. A set of coherent states is there- 
fore labeled by a quartet of complex numbers at each 
site R: z a (R) obeying the normalization conditions 
Yl a \ z a(R)\ 2 = 1- We follow the convention to indicate 
the amplitude of the | f, T> state by z\, of | |, J.> by z 2 , 
of | I, |> by z 3 and of | J,, |> by z 4 . Using this represen- 
tation, the partition function is written as 



Seff = I dt 



J2 z* a (R,t)z a (R,t)+ < z(t)\H\z(t) > 



Z= fvz a (R,t)e- s '" (20) 

in close analogy with Eq. (0). Summation over repeated 
labels a is understood. Again, by tracing out sublat- 
tice B, we reduce to an effective action defined only on 
sublattice A, formally given by: 

Seff = dt^2 z* a (R,t)z a (R,t) + J2 F i K (R,t)} 

J T^r- A TDr- D 



Re A 



ReB 



(21) 



The functional F[K(R,t)] is defined by the single site 
problem in an external field: 



e -^[K(K,t)] =Tr[//3 

f —K(R,t)U t 



(22) 



where K.(R,t) is a 4 x 4 matrix whose components K^ v 
depend on the classical field z a (R',t) defined at nearest 
neighbor sites. The explicit form of the matrix K.(R,i) 
depends on the couplings of the hamiltonian, and for the 
case we are examining is given by 

K^(R,t) = -jJ2Q(R + S,t)( l/ (R + 5,t) (23) 



where, according to our conventions, the field £ Q is sim- 
ply related to the semiclassical variable z a by: Ci = z 4, 
(2 = — 23, C3 = — z i an d C4 = z i- To lowest order in the 
spatial derivatives, i.e. taking z a (R,t) ~ Z a (R + 5, t), 
the instantaneous ground state of K.(R,t) is non de- 
generate and has components C^(R,t) with eigenvalue 
eq = —2D J where D is the spatial dimension and we 
specialized to hypercubic lattices. The three other eigen- 
values of K(_R, t) vanish and the corresponding eigenvec- 
tors are any three four dimensional vectors orthogonal to 
the ground state. Therefore the procedure outlined in 
Section II can be straightforwardly applied and requires 
the evaluation of e to second order in the lattice spacing 
a (i.e. to second order in the spatial derivatives) which 
can be obtained by standard second order perturbation 
theory: 



eo 



-2DJ + a 2 JV(z*z,)-V(z^zt) 



(24) 



Moreover the coefficient Too must be evaluated to first or- 
der in a. The lowest order just cancels the Wess-Zumino 
term, while the contribution linear in a is non vanishing 
(for smooth configurations) only in D = 1, where it gives 
rise to a residual topological term: 



Too = z a z* a + ad x {z a z* a ) + 0(a 2 ) 



(25) 



Finally, the terms Tqj (t) may be evaluated to lowest order 
in a. For slowly varying fields z a (t), only the t' ~ t region 
does contribute to the integral in Eq. (||) leading to: 






1 



2D J 



\Z a z oc — z^z^z v z v ^ 



(26) 

Combining the results ( p4]]25|]26| ) we obtain the required 
long wavelength limit of the effective semiclassical action 
for the Valence Bond SU(4) model: 



1 /■ r D+ 1 

Seff = - / dR / dt Y, E W*"! 2 + ^Q 



(27) 



The label i runs over the D + 1 space-time coordinates 
and the continuum limit has been taken in the D spatial 
directions while the imaginary time variable has been 
suitably rescaled. The coupling constant g is explicitly 
given by 



and the topological charge Q is present only in D 
where it reads: 

f 



Q = 7T~ / dxdbd x (z a z* 

Zni 



(29) 



g = a 



,D-X. 



(28) 



Tracing out a sublattice is an efficient way to take into 
account the short range antiferromagnetic correlations 
present in the model leading to an effective action de- 
scribing the much smoother fluctuations on a single sub- 
lattice. In order to support this interpretation of the 
procedure we have adopted, Fig.(@) shows the magnetic 
structure factor obtained by Lanczos diagonalization on 
a 16-site lattice: The sharp peak at wavevector k = tt 
confirms that the most relevant correlations have indeed 
periodicity of two lattice spacing. 

The formal construction of the effective action for the 
Valence Bond SU(4) model shows that the long wave- 
length and low energy physics is described by a D + 1 
dimensional NLctM or CP 3 model (with topological an- 
gle 9 = -k in D = 1). The known exact solution of 
the lattice model in one dimension implies that the CP 3 
model at g — 2 has correlations exponentially decaying 
in space-time. Moreover, the spin model is known to 
develop dimer order in the thermodynamic limit, which 
implies breaking of translational invariance by a lattice 
spacing together with breaking of parity. In going to 
the continuum limit, one sublattice has been traced out 
and then the CP 3 model remains translationally invari- 
ant but parity breaking should still occur. This picture is 
supported by Lanczos diagonalizations of the spin model 
showing that the quantum numbers of the two SU(4) 
singlet states which collapse in the thermodynamic limit 
correspond to momenta P — and P = n and opposite 
parity (i.e. reflections through a lattice site). Clearly 
the mapping we have developed neglects cut-off effects 
and holds only for sufficiently smooth configurations of 
the classical field. Therefore, the resulting estimate of 
the bare coupling constant g — 2 should be taken with 
caution but we are confident that the qualitative behav- 
ior of the CP 3 model does indeed capture the physics of 
this lattice spin system, analogously to the familiar SU(2) 
case. 

The Valence Bond model we have considered belongs 
to the SU(n)xlass already studied by AffleckEa and Read 
and SachdevcJ by use of 1/n expansion. These analy- 
sis show that the ground state breaks parity and trans- 
lational symmetry and may be described by a Valence 
Bond Solid (VBS) in one dimension, at least for suf- 
ficiently large n. The exact solution of the Valence 
Bond modelE3 confirms that this picture holds down to 
n=4. The two dimensional case is more difficult: in the 
n^ 00 limit the system has infinite degeneracy and can 
be represented as an arbitrary covering of the lattice 
by nearest neighbor valence bonds. This degeneracy is 
lifted at leading order in 1/n giving rise to a (plaque- 
tte) resonating valence bond solidEZI. In fact, the model 
maps onto a dimer hamiltoniaB-.which has been stud- 
ied by Monte Carlo techmquescil. However, diagonal- 



izations and Quantum Monte Carlo .simulations directly 
performed on the hamiltonian (|lJj)E3 suggest that the 
ground state might be a spin liquid at n=4 and therefore 
argue in favor of a phase transition between magnetically 
disordered phases as a function of n. 
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FIG. 5. Magnetic structure factor for the two SU(4) mod- 
els considered here. Lanczos data for a 16-site chain with 
periodic boundary conditions. Full dots: K = — 4J model. 
Open dots: K = 4J model. 



B. Sutherland model (K = +AJ) 

The derivation of the NLcxM appropriate for the sec- 
ond SU(4) point of the spin hamiltonian ( [l9| ) can be car- 
ried out following the same procedure adopted before and 
leading to Eqs. (BfJ) and pi]). However, the associated 
quantum mechanical problem (E21) is now defined by the 
slightly different 4x4 matrix 



K^(R,t) 



'E* 



(R + S,t)zt(R + S,t) 



(30) 



instead of Eq. (E3[). To leading order in spatial fluc- 



tuations, we can set z ll (R + 8) 



i(R). fn this case, 



the (instantaneous) ground state of the one site problem 
defined by the matrix (pO) is threefold degenerate. As 
a result, we cannot carry on straightforwardly the trace 
over one sublattice, suggesting that this procedure is not 
able to eliminate the long range oscillations in spin cor- 
relations. In turn, this means that the continuum limit 
cannot be taken just by considering the spin configura- 
tions on a single sublattice, as in usual antiferromagnets 
and larger primitive cells must be taken into account. A 
confirmation of such an interpretation comes from Lanc- 
zos diagonalizations on this model: as shown in Fig. pi 
in one dimension the spin correlations displays oscilla- 
tions characterized by the wavevector k = tt/2, implying 



a four site periodicity. This suggests a generalization of 
the procedure sketched in Section II: instead of tracing 
out one sublattice, we now keep one site every four sites of 
the chain. Therefore we need the solution of a three site 
problem, defined by the hamiltonian (fly), with "time" 
dependent boundary conditions defined by the classical 
field z a (R,t) on the two adjacent sites. Equation ( |2l| ) is 
basically unchanged, but now the sublattice A includes 
only one fourth of the sites of the chain and the free en- 
ergy functional F is defined by Eq. (g2|) in terms of a 
matrix K acting on a Hilbert space of dimension 64. To 
lowest order in spatial fluctuations of the classical field 
z a (R, t) , the full spectrum of K can be explicitly obtained 
by solving the eigenvalue equation: 

(31) 

In particular, the ground state wavefunction of the quan- 
tum problem is now non-degenerate and reads 



-^/i,^,A 



1 * 

7E Zat 



(32) 



where e aM „A is the fully antisymmetric Levi-Civita tensor. 
The corresponding eigenvalue is simply e = — 2J, while 
the Berry phase contribution gives 

Too = T z a e Q^i/A2 : /3 e /3A"^A = z a Z a (33) 

and cancels the Wess-Zumino term in the effective ac- 
tion. In order to obtain a non-trivial theory we have to 
include long wavelength fluctuations of the classical field 
z(R,t). This can be done by use of perturbation theory 
in the associated quantum three site problem. The so- 
lution is sketched in Appendix where the required terms 
are explicitly evaluated. Here we just quote the final re- 
sult, after having taken the continuum limit: The spin 
fluctuations on the sublattice are described by a NLctM 
which coincides with the one obtained for the Valence 
Bond model (p7{) at a bare coupling 



9 



6 

V5 



(34) 



larger than the estimate g — 2 obtained in one dimension 
for the Valence Bond model (J28[). Also for this hamilto- 
nian the topological angle 9 is given by 9 = it. 

As previously stressed, the Sutherland model is criti- 
cal in one dimension, with power law correlations char- 
acterized by oscillations of period 4a corresponding to a 
characteristic wavevector k — tt/2, as also confirmed by 
the Lanczos diagonalization results shown in Fig. 0. If 
the mapping onto a CP 3 theory (with topological angle 
9 — it) faithfully describes the long wavelength physics 
of the lattice model, we conclude that the CP 3 action in 
1 + 1 dimension has correlations 

< 2*(0, 0)^(0, O)z* (x, t)z a (x, t) >~ {x 2 + £ 2 )- 3/4 (35) 



On the other hand, the exact solution of the Valence 
Bond hamiltonian together with the semiclassical map- 
ping of Section IV-A implies that the CP 3 model has ex- 
ponentially decaying correlations at g — 2. These two re- 
sults may be reconciled if we assume that the CP 3 theory 
in 1 + 1 dimension undergoes a second order phase tran- 
sition as a function of the coupling g and the Sutherland 
model describes the physics at the critical point. Another 
possibility would be the occurrence of a gapless phase in 
the model, like in the celebrated CP 1 case. This alterna- 
tive explanation, however, conflicts with available analyt- 
ical and numerical evidence, pointing towards a massive 
regime at strong couplings. A phase transition in the 
CP 3 model cannot be related to a spontaneous break- 
ing of the continuous SU(4) symmetry while a possibility 
is the occurrence of parity breaking in one of the two 
phases. Of course the estimate (|34| ) we have obtained 
for the critical coupling will be renormalized by finite 
cut-off effects but the overall picture of the phase dia- 
gram emerging from the semiclassical analysis should be 
robust. This result may have implications in th* frame- 
work of the strong CP problem in field theory.c3. 



V. CONCLUSIONS 



Recent numerical simulations argued in favor of a disor- 
dered ground state in such a model for D = 2, while in 
one dimension the exact solution of the lattice hamilto- 
nian proved that spontaneous dimerization occurs. In the 
other SU(4) model we have examined, a smooth contin- 
uum limit requires a coarse graining over several lattice 
sites, suggesting that the relevant fluctuations are char- 
acterized by a wavevector k different from the antiferro- 
magnetic one. We explicitly performed the analysis only 
for the one dimensional model, where we again obtained 
the same CP 3 model which now describes fluctuations 
about k = 7r/2. The phase diagram of such a NLcrM has 
been extensively studied, particularly in I + 1 dimension, 
in connection to the CP problem in field theory. Sev- 
eral proposals have been put forward in the literature, 
including spontaneous parity breaking and deconfining 
transition. The long wavelength mapping between spin 
chains and two dimensional field theories may provide a 
clue for the final understanding of the phase diagram of 
the CP™ model, analogously to what has been found for 
the Wess Zumino Witten model. 
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In this paper we discussed some application of the 
semiclassical approach to spin systems in low dimension- 
ality. This technique is known to capture the qualita- 
tive features of quantum models and to provide a useful 
framework for the interpretation of experimental and nu- 
merical data. 

For Heisenberg chains in a staggered field we pointed 
out some difficulty and inconsistency of usual treatments, 
especially in the weak field limit, where the structure of 
the effective NLcrM turns out to be richer than expected. 
The general form of the excitation spectrum predicted 
by semiclassical approaches has been confirmed by use of 
Lanczos diagonalizations in finite clusters which also pro- 
vide some support to the usually adopted SMA. Our re- 
sults are fully compatible with existing DMRG data and 
show that different physical regimes occur in the phase 
diagram of this model. A recent DMRG investigationj-of 
aS = 2 spin chain also pointed out a similar behaviorEJ. 
The new NLcrM we derived is expected to represent only 
the low field region while, for moderate/high staggered 
fields, simple perturbative approaches, like spin wave the- 
ory, are fully adequate to describe the excitation spec- 
trum of the model. It seems unlikely that a single ef- 
fective action in the continuum limit might be able to 
encompass these two very different physical behaviors. 

A semiclassical analysis of spin-orbital models charac- 
terized by a SU(4) symmetry has also been performed 
in two different regimes. In one case we straightfor- 
wardly obtained, in the long wavelength limit, a map- 
ping to the CP 3 NLcrM describing the fluctuations of 
spin-orbitals degrees of freedom on the same sublattice. 



VI. APPENDIX 

In this appendix we briefly illustrate the procedure 
adopted for the analytical solution of the generalized 
three site problem 

(36) 

up to second order in the lattice spacing a. Here z a — 
z a + 4ad x z a + &a 2 d x z a + 0(a 3 ). We employ standard 
second order perturbation theory which gives 



Ae =< u |AK|u > 



y, \<u \AK\u n > | 2 

n#0 



(o - e, 



The unperturbed eigenvectors \u n > are the solutions of 
the eigenvalue equation (|3l]) and e„ are the corresponding 
eigenvalues. The ground state \uo > is explicitly given 
in Eq.(||) and e = -2 J. Due to the manifest SU(4) 
invariance of the eigenvalue equation, the external field 
z a may be chosen to point in the "1" direction without 
loss of generality: z a = 8 a ,i- The lowest order term gives: 

< u |AK|u >= —a 2 J [\(z* a d x z a )\ 2 - d x z a d x z* a ] 

(38) 

while the sum over excited states which appears at second 
order requires the evaluation of the matrix element 

4a J 

< U |AK|u„ >= j^^laiwUnvxdxZa (39) 
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for a generic excited state u^„i. The only solutions of 
the unperturbed eigenvalue equation (pl| ) which give a 
non vanishing contribution are those corresponding to 
the eigenvalue e — ±\/2 J and to e = 0. The former 
states are given by 1*123 = ""231 = -U132 = -U321 = 
±U3i2/v / 2 = TU2i3/"\/2 = l/\/8. The latter states are 
M123 = "321 = — U132 = — U231 = 1/2 Both excited states 
are three times degenerate: the other states being ob- 
tained by cyclic permutations of the labels (234) . Insert- 
ing these results into the perturbative expansion ( pT\ j we 
get: 



Ae = --J a 2 [\(z* a d x z a )\ 2 - 8 x z a d x z* a ] (40) 

The last required step is the evaluation of the Berry phase 
term Too to first order in the lattice spacing. Again, by 
use of perturbation theory we get 

t- / i- ^,o-T V^ < u o|AK|u„ > 

loo =< u \u > +2ilm > < u n \u > 



n^O 



U) 



(41) 



The intermediate states which contribute to the sum are 
those corresponding to e = ±V2 J which give: 



00 



z a z* a - 2ad x (z^z a ) 



(42) 



This, in turn, gives the well known topological term 
quoted in the text (|29|). 
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